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SU>CMARY 

Computational solutions have been oDtained for chemically reacting ields 

over the entire windward surface of the Space Shuttle at high 

recently devel^e^l computational method for the Space Snuttle is capable of treatl g 

three-dimensional viscous nonequilibrium air flow as well as equilibrium 

feet gas flows. A general nonorthogonal computational grid system is use 

the nLaxlsymmetric geometry. Boundary conditions take into account no^atalytic 

wall, equilibrium catalytic wall, and shock and wall slip conditions^ 

librium solutions with noncatalytic wall condition are compared to 

lytic wall solutions, the equilibrium air solutions, the perfect gas solutions 

also the Shuttle flight heating and pressure data. The comparisons show go g 

aents and correlations in most cases. 


INTRODUCTION 


Recently the nonequilibrium effect on the Shuttle reentry flowfield has oeen 
widely Investigated to reduce the surface heating by employing a proper sur ace 
aaterLl. The purpose of the present paper is to accurately predict the three 
dimensional nonequilibrium flowfield ove-r the entire Shuttle windward surrace and 
compare the result with the flight data of heating rate and pressure. a few 

typical reentry flight conditions, the nonequilibrium solutions were _ 

bo?h noncatalytic and fully catalytic wall conditions, and compared ^orreo .. 

ponding equilibrium and perfect gas solutions. The present numerical scheme was also 
extended to Include the capability to treat the nonequilibrium wall and shock slip 
conditions . 

The present numerical method (SHTNEQ) has been developed based on the two- 
•^imensionL nonequilibrium flowfield code by Miner and Lewis (ref. 1) ana the three 
dimensional perfect gas code by Szema and Lewis (ref. 2). The complete governing 
equations and the delcription of the present method are given by Kim, Swaminathan and . 
Lewis (ref 3). The SHTNEQ method uses a general nonorthogonal computational g 
sys^L to «eat the nonaxlsymm.etrlc Shuttle geometry. Since the three-dimensional 
viscous shock-laver equations are parabolic in both the streamwise and crossflow di- 
leclions, the equations are solved by a highly efficient finite-difference scheme 
developed by Murray and Lewis (ref. 4>, which requires much less computing cure than , 
?NS or^time-dependent methods. The present method can solve both subsonic and super . 
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sonic flows and requires the shock shape as initial input data. The shock shapes for 
the present Shuttle calculations were provided by the inviscid HALTS method of 
Weilmuenster and Hamilton (ref. 5). 

In the later sections, a description of the governing equations and boundary 
conditions is given, and the thermodynamic properties and chemical reaction model 
used in the present calculations are also^described . It is known that the nonequi- 
librium real gas effects persist through a wide range of the Shuttle reentry trajec- 
tory (altitudes of 122 to 50 km). In the present work, three points along the tra- 
jectory of the second Space Shuttle flight (STS-2) are chosen, and the numerical so- 
lutions are obtained over the entire windward surface of the body. These freestream 
conditions arc the same as used in ref. 3, and the present paper is an ei'tension of 
tr.e previous paper by including more calculations with comparisons to additional 
flight data and also the wall and shock slip effects. The computational resTults of 
the surface heating rate and pressure predictions are compared with the STS-2 flight 
data. The variations of some shock-layer profiles along the body are also presented. 


SYMBOLS 


C 

P 

TLXCAT 

-^1 

^2 

^3 

L 

M 

NONCAT 

NONEQL 

NSH 

P 

PG 

PHI 

py/piNF 

QW 

Re 

n 


specific heat at constant pressure 

nonequilibrium solution with fully catalytic wall 

vector in streamwlse (Cj^) direction 

vector in normal (^^) direction 

vector in circumferential (5^) direction 

Shuttle body total length, 32.84 m 

freestream Mach number 

nonequilibrium solution with noncatalytic wall 

nonequilibrium flow calculation 

shock-layer thickness nondimenslonalized by 
2 

pressure, 

perfect gas solution 

same as <\> in cylindrical coordinates 

same as p /p 
w 

. 2 

surface heating rate due to conduction and diffusion (MW/m ) 
freestream unit Reynolds number, m ^ 

dimensional Shuttle nose radius, 62.23 cm (24.5 inches) 


806 



SHTNEQ 

STS 

S/RN 

T 

f 

ref 

c 

U 

□o 

u, v,w 
U/UINF 


w . 
1 


Y/RN 
z , r , 4 ^ 
Z/L 
a 


Shuttle nonequilibrium, the present numerical method 
space transportation system 

surface distance along bodv nondimensionalized by R 

^ n 

temperature, 

2 

dimensional referer.ee temperature, U./C 

^ » poo 

Shuttle entry time from 122 km altitude 
dimensional freescream velocity 

streamwise, normal and circumferential velocity tensor components 
nondimensiona lized bv U 

CO 

streamwise velocicv, u*/U 

' ’ oo 

species production term 

body-normal distance nondimensionalized by R 

n 

reference cylindrical coordinates 

axial distance along body, same as z/L 

angle of attack, degree 

Reynolds number parameter, [u _/(c U R 

ref »oon^J 

viscosity, u*/p 

Trei 


ref 

^1 * ^2 ^^3 


reference viscosity evaluated at T 
computational coordinates 
density, p*/p 


ref 


Superscript 

* dimensional quantity 


Subscript 

i species i 

w wall value 

^ dimensional freestream value 






u 

V 


^3 


I' 
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ANALYSIS 


Governing Equations 

The governing equations are derived from the steady Navier-Stokes eouations for 
a reacting gas mixture as given by Bird et al. (ref. 6) , and they are written in a 
surface-oriented general nonor thogonal coordinate system (see fig. 1). The Ci co- 
ordinate consists of straight lines in the surface-normal direction. At the body 
surface, the coordinate is chosen to coincide with the coordinate of the ref- 
erence cylindrical coordinate system. The coordinate system requires orthogonalitv 
on^y at the cady surface. The normal velocity v and normal coordinate ^2 assumed 
to be the order of e, and all terms which are of higher order than r. are neglacted ii 
the governing equations. The local physical velocity vector is defined as 

where the u, v and w are tensor velocity components in the computational coordinate 
system. Only laminar flow is considered in the present analysis. The derived nondi- 
raensional form of the three-dimensional viscous shock-layer equations for a reacting 
gas mixture is given in ref, 3. 


Boundary Conditions 

At the body surface, the slip and temperature- jump boundary conditioHw can be 
used if necessary. The nonequilibrium wall slip equations given by Hendricks (ref. 7) 
have been rewritten for the present coordinate system and chemical oodel. The cal- 
culated Reynolds number parameter e was less than 0.108 for the present test cases 
which indicates that the slip effects on the heating rate and surface pressure will 
be small (see e.g. ref 8), but the slip effects on some shock-layer profiles over 
the nose region are calculated and presented. The wall temperature is specified by 
the STS— 2 flight thermocouple data. In the present calculations, the wall species 
concentration is dictated by the noncatalytic or fully catalytic condition, but the 
boundary condition can easily be extended to include the effects of finite wall 
catalycity on the recombination of dissociated air. At the low surface temperature 
of Shuttle, the equilibrium catalytic wall condition can be replaced by the fully 
catalytic wall condition. 

In the present method, shock shape information is necessary as an input which is 
used for the calculation of the shock— boundary condition. The three— dimersional 
shock-boundary conditions with slip effects (modified Ranklne— Hugoniot jump relations) 
given by Murray and Lewis (ref. 4) have been extended to include finite— rate chemistry 
and the nonorthogonal coordinate system. Two-dimensional shock-normal coordinates 
are defined in the plane which contains both the freestream velocity vector and the 
vector which is normal to the local shock surface- Then, the freestream velocity 
vector is written in the shock-normal coordinates, and two-dimensional shock-crossing 
conditions are calculated in the shock-normal coordinates. The known after-shock 
quantities are rotated into the three-dimensional computational coordinate^. 


Thermodynamic and Transport Properties 

Multi-component ionizing air is considered to be a mixture of then., lly perfect 
gases, and the thermodynamic and transport properties for each species are calculated 
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using the local temperature. The properties for the gas mixture are then determined 
in terms of the individual species properties. The enthalpy and specz^zz^c heat of 
each species are obtained from the thermodynamic data tabulated by Brcvne (refs. 9-11). 
A second-order Lagrangian method is used to interpolate the values at = given tempera- 
ture, The viscosity of the individual species is calculated from the cnrve fit re- 
lation given by Blottner (ref. 12), and the thermal conductivity of each species is 
calculated from the Eucken semi-empirical formula using the species vrscosity and 
specific heat. After Che viscosity and thermal conductivity of the inza.vldual species 
are calculated, the viscosity and thermal conductivity of the gas mixt^ra are calcu- 
lated by the method suggested by Armaly and Sutton (refs. 13, 14). In cae present 
work, the diffusion model is limited to binary diffusion with the binary diffusion 
coefficient specified by the Lewis number of 1.4. 


Chemical Reaction Model 


It is assumed that the reactions proceed at a finite rata, and the rate 

of production terms wi of the individual species are included in the energy equation 
and the species continuity equations. The Wj^ terms are functions of born the tem- 
perature and the species concentrations, ard they must be rewritten so rnac the tem- 
perature or the species concentrations appear as one of the unknovms as given in 
ref. 1. In the present calculations, the chemical reaction model and rbe reaction- 
rate constants are taken from Blottner (ref. 15). Seven (7) chemical species are con- 
sidered in the reactions; viz., 0, O 2 , ^0, N, NO*^, N 2 , and e , The foLlewing pure air 
chemical reactions are used for the present study: 


1. 

02 


Ml 

->• 

A- 

20 

+ 

Ml 

2. 

N2 


M2 

■*- 

2N 

+ 

M2 

3. 

»2 


N 

Ar- 

2N 


N 

4. 

NO 

+ 

M3 


N 

+ 

0 - 

5. 

NO 

+ 

0 

-> 

O 2 

+ 

N 

6. 

N2 

+ 

0 

-► 

A- 

NO 


N 

7, 

N 

4- 

0 


NO*-+ 



where Ml, M2 and M3 are the catalytic third todies (i 
duction terms are for nonequilibrium flows, the prese. 
in obtaining a converged solution whenever the flow c 
The difficulty is severe, particularly at the stagnate 


Since tae rate of pro- 
od encountexs difficulty 
ns approach equilibrium, 
point . 


Numerical Soltition 

Davis (ref. 8) presented an implicit finite-difference method to solve the vis- 
cous shock-layer equations for axially symmetric flows, and Murray and l^ewis (ref. 4) 
further developed the scheme for three-dimensional flows. In the present work, the^ 
method is extended to the chemically reacting three-dimensional flcwfield solution in 
a surface-oriented nonorthogonal coordinate system. Since the viscous shock-layer 
equations are parabolic in both the stxeamwise and crossflow directions, the equations 
are solved by a highly efficient finite-difference scheme- The continiiiry and normal 
momentum equations are solved in a coupled form to promote convergence- The shoe* 

stand-off distance is evaluated by integrating the continuity equation. 


The solution begins on the spherically blunted nose by obtaining an axi symmetric 
solution in the wind-fixed coordinate system. The axisymmetric solution is rotated 
into the body-fixed coordinates and is used as the initial profile for rhe three- 
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ditnensional solution. The three-dimensional solution begins in the windv:-.rd 
and marches around the body obtaining a converged solution at each ^^-stec". After 
completing a sweep at a Cl”niarching station, the procedure then steps dowTr^tream in 
E^l and begins the next 53 -sveep. At each point the equations are solved An the fol- 
lowing order: (i) species continuity, (li) g^-momentum, (iii) energy, ( i*r) 

turn, (v) integration of continuity for shock-layer thickness, and (vi) cc*uplec ern- 
tinuity and normal momentum equations. 


RESULTS -AM) DISCUSSION 

In order to predict the Shuttle reentry flowfieid, three test cases -^ere znrsen 
and the viscous windward flowfieid solutions were obtained using various tnentza- 
models. For the wide range of the Shuttle reentry conditions (above 50 knr alttt^ey 
the nonequil ibr iuii effects can occur, and the nonequilibrium effects are -.argest near 
the nose of the bedy and around t = 450 sec on the trajectory (ref. 16). The a_tt- 
tudes selected for the presen^ calculations are 81, 70 and 60 km (t ~ 250. 460, tjO 
sec, respectively). Detailed freestrean conditions for the three test cashes are 
given in table I. The inviscid input shock ^shapes for the present calculanions have 
been provided by the HALIS method for an angle-of-attack of 40 deg for boAh perrett 
gas and equilibriuia air. The inviscid HALIS shock was available only up zs> z/1 = 

0.5 or less. The viscous flowfieid solutions for perfect gas and equilibirluin air 
been obtained up to z/L = 0.5 in order to compare with the nonequi librium solunlrns. 

The nonequil ibriun solutions, however, were obtained for the entire windwacrd surrace 
up to the body end using an extended shock. The shock extension was done using tae 
STEIN (ref. 17) solution of the shock shape for an angle-of-atcack of 25 Cteg. Tne ex- 
tended shock was scaled and smoothed before being used as input data. In order to 
enhance the accuracy of the nonequilibriinn viscous solution, a global iterraticn nas 
been performed using the viscous output shock as an input. 

The nonequilibrium solutions have been obtained for both noncatalytic and lu.ly 
catalytic wall conditions for the purpose of comparison. The cross-sections 01 toe 
modified Shuttle orbiter which have been used for the inviscid and viscous solntims 
are depicted in fig. 2. Presented results include the surface heating fate, surtace 
pressure, shock shapes, a few shock-layer profiles, and finally slip effeetts over the 
nose region. The flight heating— rate data obtained from the method by Throckmertm 
(ref- 18) are used for comparison ^/Ith the computational results. The fli^t -reasurer— 
ment data "of pressure are used for comparison with the present surface-pressure pre- 
dictions. 


Surface Heating Rate 

The heating-rate predictions along the windward centerline are compareid with 
each other in fig- 3 for Case 1. The nonequilibrium solution with noncata_ytic vail 
condition agrees well with the flight da^a for most of the region. The ncmequi:.itriunr 
solution with fully catalytic wall condition shows quite close agreement wlxh tne 
equilibrium =^ir solution. The heating-rata prediction from the perfect gas mocei is 
below the equilibrium air solution but well above the noncatalytic wall soluticn for 
the entire body. The reason for the local mismatch arcjuno t:/L — C.c is noc currently 

known. At z/L = 0.4, the perfect 'gas solution is 20% lower chan the equil iorlur* so- 

lution and 50% higher than the noncatalytic wa7 1 solution. In fact, the s^rzace 
finite catalytic effect for Case 1 is negligible compared to that for CaS'es 2 and 0 , 
Gue to Che altitude dependence cf Che surface catalytic activity (ref. 19). Tne de- 
crease of the heating after z/L = 0.8 is due to the slope change of the beefy suriaze* 
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For Cases 2 and 3, the g»;neral trends of the computed heating rates are similar 
to the result of Case 1, but the noncatalytic wall solution underpredicts comparer to 
the flight data especially on the nose region and the body-end region as shown it 
figs. 4 and 5. The discrepancy over those regions may be due to the surface cata^ytii 
effect as recently discussed by other investigators (refs. 16, 19 and 20). In - 

the nonequilibr iura solution with noncatalytic wall condition along the body is stcvn 
for various 4^-planes, together with the corresponding equilibrium air solution r_r 
Case 2. The equilibrium air solution is much higher than the ncncatalytic wall seda- 
tion for all the i-planes. In fig. 6, the surface heating distributions along tne 
body for £tll the z— planes are shown. The sudden increase of the heating rate at 
(() = 80 deg and z/L = 0.5 is due to the spanwise slope change of the body surface 
along the oody (see fig. 2). At z/L - 0.6 and - 90 deg, the solution did not cci- 

verge, due to the severe surface slope change of the wing tip section- The calcu- 
lated windward spanwise heating rates at two axial stations (z/L = 0.2 and z/L = 
for Case 2 are shewn in fig. 7, The comparisons among the various chemical models 
show similar trends and correlai: Lons for the spanwise heating rate distributions cue 
to the flow expansion around the body. The flight data were available only at tie 
windward centerline (i> = 0 deg), and the data agree well with the noncatalytic va^_ 
solution . 


Surface Pressure Comparison 

In fig. 8 the surface-pressure distribution over the entire Shuttle windward 
surface is presented together with the available flight data for comparison. The 
cj>-planes from 10 deg to 40 deg were omitted on the plot because the results for those 
planes were almost identical with the result for the windward centerline. The agree- 
ment with flight data is good especially on the windward centerline for Cases 2 arc 3- 
In fig. 9 the spanwise surface pres-sufe distributions are shown at two axial staitcas 
(z/L - 0.1 and z/L = 0.2)./ The present calculation tends to underpredict over cne 
flow expansion region (off the centerline) compared to the flight data (e.g. about 2G.« 
underprediction at b ~ 67 deg, z/L = 0.1). This disagreement may he cue to the un- 
certainty of the flight data. A computation using a stialler (J)-stepsize may also re- 
duce the discrepancy. The present solution by the SHTXEQ method used b-stepsize ci 
10 deg around the body, and this stepsize may not be small enough for the nonclrcr^ar 
cross section of the Shuttle geometry (see fig. 2). If more 45-pIanes than the presenr 
10 planes were included in the computation, the current relatively large storage re- 
quirement and computing .ime would increase accordingly. 


Shock-Layer Thickness 

A comparison of the shock-layer thickness distributions along the body at the 
{p = 0 plane which has been obtained from various chemical models is shown in fig. -0. 
All the viscous shdek shape results except the inviscid shock are ^rom the once 
globally iterated results. Wnen an inviscid input shock is not very accurate, the 
output shock shape is, in general, different from the input shock- In such a case, 
a global iteration is necessary in order to refine the entire flawfield solutioa. 
Thus, for the present three test cases, all the presented viscous flo'wfield soluritns 
are from the first global iteration. The inviscid shock is from the H.ALIS code anc 
was available only up to z/L = 0.5 as mentioned earlier. The inviscid HALTS shock 
for OL = 40 deg has been extended by the inviscid shock of ct = 25 deg irora the 
method. The extended shock was scaled and smoothed, and then used as the input stock, 
data for the initial calculation (zeroth iteration). When the various viscous shec-k- 
layei thicknesecs ar^^ compared to the inviscid perfect gas shock at z/L = 0.4, the 
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viscous perfect gas shock is 86.7%, the nonequilibrium shock is 52.8% and the 
ecuilibrium shock is 36.1%. The shock from the fully catalytic wall condition — s 
almost identical to the noncitalytic wall Solution- 


Comparison of Shock“Layer Profiles 

The nonequilibrium flowfield structure of the viscous shock-layer at a few 
selected axial stations on the windward centerline is depicted for both the non- 
catalytic and fully catalytic wall conditiuas in figs. 11-14. ihe profiles inc-ice 
temperature, tangential velocity and mass fractions of oxygen and rnTrogen atoms icr 
Case 2 (t = 460 sec). The temperature and velocity prcfilcs in figs. 11 and 12 rhcw 
that the viscous effects are dominant across the entire shock-layer, especially 
the forward part of the body. At the axial station of z/L = 0.046 j-ti the plot -- 
temperature profile, the wall temperature gradient of the fully catalytic w-all s:iu- 
tion is larger than that of the\oncatalytic wall solution, which produced a 34^ 
larger conduction heating rate. At the same station the heating rate due to 
mass diffusion was 93% of the conduction heating rate in the fully catalytic wal^ sc— 
lution. The diffusive heating in the noncatalytic wall case was, of course, negligiht 
Tae surface chemical catalycity has negligible effect on the velocity profile as 
shown in fig. 12. The mass-fraction profiles of oxygen and nitrogen clearly show 
effects of the noncatalytic and fully catalytic wall corditions. In the ncncata^jti-r 
vail case, the oxygen atom concentration at the wall remains almost constant alcrrs 
the body, while the nitrogen atom concentration is reduced downstream due tc mori 
rap id recomb ina tion. 


Slip Effects 

The SHTNEQ method has been further extended to include the shock slip and tie 
vail slip condition's. In high altitude freestream conditions, the conventional 
frozen shock crossing of Rankine-Hugoniot relations for nonequilibrium flows givrs 
poor prediction of the after-shock quantities. It is known chat the slip effects oit 
surface-measurable quantities like heating rate and pressure are significant, 
especially for reentry bodies with a small nose radius at very high altitude. Frr 
the Space Shuttle geometry at the test case freestream conditions, however, the xil- 
culated Reynolds number parameter c was less than 0.108 which indicates that the slxp" 
effects on the surface-measurable quantities will not be significant (ref. S) . In 
fig. 15, the shock-slip effects on the temperature profile and mass fraction of 
oxygen and nitrogen atoms are shown at the stagnation point. The slxp Lemperature ac 
the shock is less than the no— slip temperature by 1500 K. The shock— slip effect 
oa the oxygen mass-fraction distribution across Che shock-layer is limited to the 
region near the bow shock, but foi nitrogen the shock-slip effect is propagated c-1 
through the shock— layer. Figure 16 shows the wall— slip effects on the surface L-sn— 
perature and axial flow velocity jumps over 'the nose region for Case 2. The amomt 
of the temperature jump is about 200 K at the stagnation point, and the slip 
velocity is 0,0056 times the freestream velocity at S/RN = 0.8. 


Computing Times 

The computing times required for the flowfield computations of all the test 
cases are liszed in table II. The computing times are based on an IBM 370/3081 
general purpose computer. The nonequilibrium computations took about one and a half 
hours CPU time for solving the entire Shuttle windward surface. When an input shick 
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data is not accurate, generally a global iteration xs lequired, thus consuming mere 
computing time uhan the tabulated one. The computing times for the perfect gas and 
equilibrium air are for the solution of the first half of the body (up to z/L = O.o). 
The solution of the perfect gas or equilibrium air flows took relatively small cen— 
puting times (less than 20% of nonequilibrium case^ . Tl.e axial marching step—sizc^ 
are controlled internally in the code considering the number of local iterations 
taken. A fixed input of 51 or 101 grid points was used in the surface-normal di- 
rection and 10 planes were used around the body for the windward surrace (10 deg 
step— size) . The leeward surface the Space Shuttle was not considered, because a 
solution cannot be obtained by the present method cue to strong ilow separation. ^ The 
storage requirement of the present SHTNEQ code is x52 kilo-bites in the i.3M 370/^-31 
computer . 


CONCLUSIONS 

In general, the ron^pur \ -lal results of surface heating rate'for the three- 
dimensional nonequiUb) ^ -afield over the Space Shuttle compare well with the 

available flight data- The flight heating rate data are higher than the loncata-^ 
lytic wall solution especially on thz nose region and the body end due to the surface 
finite catalytic effect. The nonequilibrium solution with fully catalytii wall gives 
quite close agreement with the equilibrium heating rate prediction. The perfect gas 
solution of surface heating rate is less than the equilibrium solution but higher 
than the flight data for the entire region. The calculated pressure distribution also 
shows good agreement with the flight data. The calculated moneq’iilibriun shock- and 
wall-slip effects on the heating rate were negligible for the present test cases. The 
computing times taken for the nonequilibrium calculation are reasonable conslriericg 
the large size of the computational grid due to the complex Shuttle geometry and ,r.e 
chemical reactions of seven species. Further work is planned to calculate the ef.tcts 
of finite catalytic wall conditions on the three-dimensional nonequilibrium flowfieid. 
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TABLE I. TEST CASE FREESTREAM CONDITIONS 


Case 

No. 

t 

(sec) 

Altitude 

(km) 

a 

(deg) 


Reoo 

(m-l) 

UoD 

(km/ s) 

^oo 

CO 

(atm) 

1 

250 

35.74 

41.0 

26.6 

2726, 

7.53 

199. 

3.587E-C6 

2 

460 

74.98 

40.0 

25.5 

15686. 

7.20 

193. 

2.142E-03 

3 

650 

71.29 

39.4 

23.4 " 

*25756. 

6.73 

205. 

3.965E-05 



TABLE 

II. 

COMPUTING 

TIMES® 

FOR TEST CASES 




Case 

z/L 


Grid Size 

of 


CPU Time 

Flow Model 

- No. 

from - to 

-steps C--pts 

±. L 

^3- 

-planes 

(H:M:S) 

Nonequilibrium 

\ 

0. 

- 0.93 

111 

51 


10 

1:20:34 

Noncatalytic 

2 

0. 

- 0.93 

129 

51 


10 

1:50:28 

3 

0. 

- 0.93 

142 

51 


10 

1:51:40 

Non equilibrium 

1 

0. 

- 0.93 

123 

51 


10 

1:23:02 

Fully Catalytic 

2 

0. 

- 0.93 

124 

51 


10 

1:26:53 


3 

0. 

- 0.93 

124 

51 



1:27:22 


1 

0. 

- 0.50 

78 

101 

■ 

10 

0:07:07 

Perfect Gas 

2 

0. 

- 0.50 

78 


■ 


0:07:0b 


3 

0. 

- 0.50 

78 

mm 

■ 


0:07:05 


1 

0. 

- 0.453 

72 

101 


10 

0:08:32 

Equilibrium 

2 

0. 

- 0.4^53 

72 

101 


10 

0:08:28 


3 

0. 

- 0.453 

72 

101 


10 

0:08:27 


^CPU time cn IBM 370/3081, H'-OPT2 comp'Ier 






ORIGINAL FAGS S 
OF POOR QUALITY 


o 



Figure 3.- CoEparison of measureu and calculated heating 
races along the windward centerline of Case 1, 


□ 



Figure 4.- Comparison of measured and calculated heating 
rates along the windward centerline of Case 2. 
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Figure 5.- Comparison of measured and calculated beating 
rates along the windward centerline of Case 3. 




(a) 
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Figure 6-- Surface heating rate distributions along the body 
for different v”planes of Case 2. 
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(b) Case 2. 

Figure 8.- Comparison of surface pressure distributions along the body. 
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Figure 10.- Comparison of globally iterated shock laye: 
thickness solutions at (J) - 0 for Case 1. 



Figure: 


11.- Shock layer profiles of temperature at various 
axial body locations along <{) = 0 for Case 2. 
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Figure 14.- Shock layer profiles of nitrogen mass fractio 
various axial body locations along = 0 for Case 2. 


o 



Figure 15.- Shock slip effect on the profile of t-jiiiperat 
and mass fraction at stognation point of Case 2. 










